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Inelastic neutron scattering experiments are performed on a nanoporous metal complex Cu-Trans- 
1,4-Cyclohexanedicarboxylic Acid (Cu-CHD) adsorbing O 2 molecules to identify the spin model of 
the 02 -based magnet realized in the host complex. It is found that the magnetic excitations of Cu- 
CHDs adsorbing low- and high-concentration O 2 molecules are explained by different spin models, 
the former by spin dimers and the latter by spin trimers. By using the obtained parameters and also 
by assuming that the levels of the higher energy states are reduced because of the non-negligible 
spin dependence of the molecular potential of oxygen, the magnetization curves are explained in 
quantitative level. 

PACS numbers: 75.50.-y, 75.50.Xx, 78.70.Nx 


I. INTRODUCTION 

Natural oxygen, the second abundant constituent in 
the air, is a magnet having spin S=\ induced by a couple 
of 7Tg electrons. The oxygen molecules are condensed and 
crystalized at 54 K by van-der Waals interaction. The 
solid oxygen realized at the low temperatures include a 
phase with monoclinic C^jm at T < 23 K, /3 phase with 
hexagonal i?3m at 23 K < T < 43 K, and 7 phase with 
cubic Pm3n at 43 K < T < 54 K.[I|, 0 Since the mag¬ 
netic system is phenomenologically coupled to the lattice 
system thorough magneto-elastic coupling, the successive 
phase transition is magneto-structural transition accom¬ 
panied by the change of the spin Hamiltonian. Magnetic 
susceptibility measurements suggested that quasi-two- 
dimensional spin system in the a phase, triangular spin 
system in the ff phase, and one-dimensional spin system 
in the 7 phase were realized, [ij implying pressure induces 
more phases including metal, Q superconductivity, Q, 
and an insulating magnetically ordered phases. Q The 
variety of the phases suggests that the energy scales of 
the coupling among multi-degrees of freedom, i.e., lat¬ 
tice, spin, and orbital, are close to each other in the soft 
solid crystalized by Van der Waals interaction. Indeed 
the molecular potential between a pair of O2 molecules 
is strongly dependent on the spin state, which means that 
the geometrical configuration of the O2 molecules and the 
spin state are closely correlated. [iSIl To manipulate the 
oxygen molecule O2 and to artificially synthesize a novel 
type of 02-based magnet is, therefore, the challenge in 
the new field of magnetism. 

Pioneering work is found in the adsorbed O2 on the 
surface of graphite, where triangular lattice of O2 is 
realized. [1^ [1^ Combination of the magnetic suscepti¬ 
bility and neutron diffraction measurements reveals the 
magnetically ordered state at low temperatures. An¬ 
other root for the 02-based magnet is to utilize nano¬ 
materials such as microporous metal complexes. 


nanoporous silica, [2^ [2l| or carbon nanotubes. The 
adsorbed molecules staying at the minimum of the Van 
der Waals potential in the nanopore form a supercrystal, 
leading to the realization of the 02-based magnet. 

We focus our attention on a metal complex 
having one-dimensional nanopores, Cu-Trans-1,4- 
Cyclohexanedicarboxylic Acid Cu2(OOC-C6H4-COO) 
abbreviated as Cu-CHD. [I^ The adsorbed O2 molecules 
form dimer-like structure in the nanopores. [23| Indeed, 
the magnetic susceptibility of Cu-CHD adsorbing 0.18 
mole of O2 per half of formula unit showed a rapid 
decrease with the decrease of the temperature, which is 
consistent with a spin dimer model. It was, however, 
explained not by S'=l spin dimer but rather by S' = 
1/2 spin dimer. To explain the unusual bulk prOTerty, 
the spin-dependent Van der Waals potential p-Q was 
discussed, and the anomalous energy spectrum was 
proposed. [13 The precise energy scheme is to be directly 
revealed by a spectroscopic technique. 

The magnetic susceptibility of Cu-CHD adsorbing high 
concentration of O2 molecules showed qualitatively dif¬ 
ferent behavior from that of the low concentration one. 
Curie-Weiss like behavior was enhanced in the low T 
region, [l^ The result suggested that the spin model and 
the corresponding energy scheme can be tuned by the 
concentration of the O2 molecules in the 02-based mag¬ 
net. The detailed model is to be experimentally identi¬ 
fied. 

In the present study we carried out inelastic neutron 
scattering (INS) experiments at low temperatures in or¬ 
der to observe the low-energy excitations and to clarify 
the spin model of the 02-based magnet realized in Cu- 
CHD. We found that the magnetic excitation of Cu-CHD 
adsorbing less O2 is explained by a spin-dimers model of 
which the exchange constants are normally distributed. 
In contrast, that of Cu-CHD adsorbing more O2 is ex¬ 
plained by a spin-trimers model with the distributed ex¬ 
change constants. It was found that the spin model is 


2 


tuned by the concentration of O 2 in this system. Based 
on the spin models with the conclusive parameters de¬ 
termined by INS experiments and by introducing the re¬ 
duction of the energy levels of higher states, we quanti¬ 
tatively explained the magnetization curves. The effect 
of the spin-dependent Van der Waals potential to the 
anomalous energy scheme is experimentally confirmed. 


II. EXPERIMENTAL DETAILS 

Polycrystalline Cu-CHD was synthesized in methanol 
solution. The details are described elsewhere. The 
total mass was 5.93 g. Since Cu-CHD includes many 
hydrogen atoms that have large incoherent neutron scat¬ 
tering cross sections, we paid special attention to extract 
the net contribution of the adsorbed oxygen molecules. 
The first step was to prepare the bare Cu-CHD sample 
in which any type of guest molecules including H 2 O, N 2 , 
O 2 , etc., were eliminated. The second step was to mea¬ 
sure the INS spectrum of the bare Cu-CHD sample as 
background data. The third step was to prepare the oxy¬ 
gen adsorbed Cu-CHD sample. Here we used exactly the 
same sample that was used for the background measure¬ 
ment in the second step. The forth step was to measure 
the INS spectrum of the oxygen adsorbed Cu-CHD. Ide¬ 
ally the net contribution of the adsorbed oxygen can be 
extracted by subtracting the INS spectrum in the sec¬ 
ond step from that in the forth step. In reality we found 
additional intensity that increased with the increase of 
the wave number transfer Q even after the background 
subtraction process. In addition we found flat intensity 
that was independent both on Q and the energy transfer 
Hlo. During the analyses on the data after background 
subtraction, therefore, we used additional background of 
which the form was A{fujj)Q^ -|- B. Here A{huj) and B 
were free parameters. We could not identify the origin of 
the additional background, even though we could make 
some speculations; A{huj) term would be from the defor¬ 
mation of the host compound after the oxygen adsorp¬ 
tion, and B term would be from the difference of white 
background such as electronic noise between before and 
after the oxygen adsorption. 

The procedure of preparing the bare Cu-CHD is as fol¬ 
lows. The Cu-CHD sample in the Aluminum made con¬ 
tainer was warmed at about 105 °C and was evacuated by 
a vacuum pump with monitoring the pressure. We used 
the gas evacuation/introduction system described in Fig. 
2 in RefUi We waited until the pressure reached a few 
Pa, which is the capacity limit of the vacuum pump, and 
we confirmed that the Cu-CHD sample had discharged 
all the molecules. We, then, seal the container by clos¬ 
ing a VI valve in Fig. 2 in Refill, put the sample in 
a liquid helium cryostat, cool the temperature down to 
4.5 K, and collected INS spectrum of the bare Cu-CHD. 
After the background measurement we set the tempera¬ 
ture of the bare Cu-CHD at 110 K and introduced O 2 
gas in the nanopore by using a buffer container of which 


the capacities were known. The pressure of the gas was 
kept below 8x10^ Pa. The volume of the introduced 
oxygen was estimated by the capacity of the buffer con¬ 
tainer and the pressure of the O 2 gas inside the con¬ 
tainer. After the O 2 introduction procedure the sam¬ 
ple container was sealed, and was cooled down to 4.5 
K. All measurements were carried out at T=4.5 K. We 
prepared two samples that adsorbed different amount of 
O 2 in order to examine the O 2 concentration dependence 
of the spin model realized in the nanopores. Exactly 
the same polycrystalline Cu-CHD was used as the host 
compound for both samples. In prior to the preparation 
of each sample we discharged the gas in the nanopore 
at 105 °C according to the procedure described above. 
The O 2 mole ratios sealed in the sample container were 
02/(Cu atom)=0.3±0.07 and 2.0±0.4, which are labeled 
as O. 3 O 2 -(Cu-CHD) and 2.OO 2 -(Cu-CHD), respectively. 

Neutron scattering experiment was carried out using 
the cold neutron ch opp er spectrometer AMATERAS in¬ 
stalled at J-PARC.[^ [ 2 ^ The initial neutron energies 
i5i=3.135 meV and 15^=7.743 meV were used, and their 
energy resolutions are about 0.15 and 0.45 meV at the 
elastic position. 


III. RESULTS 

Figures [Ija) and (b) show the raw inelastic neutron 
scattering (INS) spectra of bare Cu-CHD and O. 3 O 2 - 
(Cu-CHD), respectively. The initial energy Ei is 7.743 
meV. Strong incoherent scattering from hydrogen atoms 
and/or functional groups were observed in the wide re¬ 
gion in the energy transfer {hu) - wave number trans¬ 
fer (Q) space in both panels. Figure [IKc) shows typical 
one-dimensional (ID) energy cut obtained by integrating 
the intensity in the range ofl.2A ^ < Q < 1.4A ^ 
for the raw spectrum of 0.3O2-(Cu-CHD) and that of 
the bare Cu-CHD. Statistically meaningful increases are 
observed at about 2.5 meV and 4 meV in 0.3O2-(Cu- 
CHD). The enhancement of the intensity is due to the 
additional scattering of the adsorbed oxygen molecules. 
The decrease of the intensity is observed at about 1 meV 
in Fig. [TJb). As we described in the Section H the INS 
spectrum of the bare Cu-CHD is regarded as background, 
but its intensity is comparable to the 0.3O2-(Cu-CHD) 
because of the large incoherent scattering from the hydro¬ 
gen atoms included. Even though the careful background 
measurement, the background has rather larger intensity 
than 0.3O2 -(Cu-CHD), of which the origin would be un¬ 
avoidable artifacts. 

Figure [IJa) shows the INS spectrum of 0.3O2-(Cu- 
CHD) after the background subtraction in the low-energy 
region for the data set of 75^=3.135 meV. The absence of 
excitation at fio; ;< 2 meV is confirmed in the experimen¬ 
tal error. Figure [2jb) shows the INS spectrum of O. 3 O 2 - 
(Cu-CHD) in the high energy region for the data set of 
£1^=7.743 meV after the background subtraction. The 
dispersionless excitations are observed at Iia;~2.4 meV 
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(a) Bare Cu-CHD raw data 
F..=11A^ me.V 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Q (A->) 


(b) 0.3O2 -(Cu-CHD) raw data 




FIG. 1: (Color online) (a) Contour map of the raw neu¬ 
tron spectrum bare-(Cu-CHD) by using iJi=7.743 meV. (b) 
Contour map of the raw neutron spectrum 0.3O2-(Cu-CHD) 
by using i5i=7.743 meV. (c) One-dimensional energy cut ob¬ 
tained by integrating the intensity in the range of 1.2 A < 
Q < 1.4 A ^ for raw data of 0.3O2-(Cu-CHD) and bare Cu- 
CHD. 


and 4 meV. These two excitations are also probed as two 
peaks in the ID-energy cut after the background subtrac¬ 
tion shown in Fig. He). 

To examine the (^"dependence of the excitation at 
tvjj~4: meV, the lD-(5 cut obtained by integrating the 
intensity in the range of 3.5 meV < huj < 4.5 meV is 
shown in Fig. H^)- The intensity of the excitation has 
a kink at about 1.2 and it decreases with increas¬ 
ing Q. The dispersionless behavior means that the origin 


of the excitation is a cluster. The suppressed intensity 
at large Q means that the excitation is dominated by a 
magnetic scattering. The excitation induced by the oxy¬ 
gen adsorption is, therefore, the magnetic one from some 
types of the O 2 molecules cluster. The kink at finite Q 
means that the spin correlation is antiferromagnetic, and 
the maximum Q exhibits the inversed length-scale of the 
intra-cluster distance. As for the excitation at fiw~2.4 
meV in Fig. Hb), the intensity shows monotonic de¬ 
crease with the increase of Q in the range of (5 < 2.2A 
and, then, it increases in the range of 2.2A < Q. We 

considered that this excitation is not intrinsically mag¬ 
netic one, which will be discussed later. 

Raw INS spectra of the bare Cu-CHD and the 2 .OO 2 - 
(Cu-CHD) with Ei = 3.135 meV are shown in Fig. [3] (a) 
and (b), respectively. A flat excitation a.i hw ^ 0.4 meV 
is newly observed in the 2.0O2-(Cu-CHD) in addition to 
smeared intensities in wide huj - Q space. ID-energy cuts 
of the 2.0O2 -(Cu-CHD) and the bare Cu-CHD in the 
integration range of 1.0 A“^ < Q < 1.2 A“^ is shown in 
Fig.[3](c). Peak structure at 0.4 meV is clearly observed. 
Raw INS spectra of the bare Cu-CHD and the 2 .OO 2 - 
(Cu-CHD) with Ei = 7.743 meV are shown in Fig. [3Dd) 
and (e), respectively. The spectrum in Fig. [3Kd) is the 
same data as shown in Fig. [TJa) but the intensity scale 
is different. Broad excitations at 3 meV ^ fiuj < 5 meV 
in low Q region are newly observed in 2.0O2-(Cu-CHD). 
The ID-energy cuts for Ei = 7.743 meV is shown in Fig. [3] 
(f). The Cu-CHD adsorbing more O 2 molecules exhibits 
the large enhancement of the neutron intensity compared 
to those of 0.3O2 -(Cu-CHD). 

The INS spectrum of 2.0O2-(Cu-CHD) by using the 
£1^=3.135 meV after the background subtraction is shown 
in Fig. The dispersionless excitation is observed at 

i?~0.4 meV. The ID-energy cut in Fig. HJe) exhibits that 
the energy width of the excitation is almost the resolution 
limited. The INS spectra of 2.0O2-(Cu-CHD) by using 
the Ei=7.7A3 meV after the background subtraction is 
shown in Fig. Hb)- The dispersionless excitations are 
observed at fia;~2.4 meV and 4 meV. These excitations 
have the broad width as shown by the ID-energy cut in 
Fig. Hd). 

In order to examine the (^-dependence of the excitation 
at 0.4 meV, the ID-Q cut obtained by the integration 
in the range of 0.3 meV < E < 0.5 meV is shown in 
Fig. He). The dispersionless excitation has broad peak 
at (5=1.2 A“^ of which the wave number corresponds 
to the distance between O 2 molecules. This excitation 
is not clearly observed in 0.3O2-(Cu-CHD). The energy 
of the excitation is quite close to the anisotropy energy 
of the natural oxygen in the gas phase reported as 0.5 
meV.jl^ This suggested that the excitation is due to the 
single-ion-anisotropy of the adsorbed oxygen molecules. 

The lD-(5 cut for the excitation at 4 meV is obtained 
by the integration range of 3.5 meV < E < 4.5meV 
in Fig. HO- The excitation has the peaks at about 
1.2 A“^ and the intensities decrease with increasing Q. 
The qualitative behavior is the same as that of O. 3 O 2 - 
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(a) 0.3O2-(Cu-CHD) subtracted data 

.135 meV 


(b) 0.3O2-(Cu-CHD) subtracted data 
E=1.143 meV 
30 


1.0 

e(A-‘) 


1.5 „2.0 

Q (A-‘) 


(c) 



(d) 



FIG. 2: (Color online) (a,b) Contour map of the neutron spectrum of 0.3O2-(Cu-CHD) after the background subtraction by 
using (a) i?i=3.135 meV and (b) Ei—7.743 meV. (c,d) Neutron intensity for 0.3O2-(Cu-CHD) after the background subtraction 
by using Ei—7.743 meV. The green cnrve is the nentron cross section of a gronp of spin dimers having distributed J. The 
blue curve is the additional background having Q dependence. The red curve is the sum of green and blue curves. The dashed 
line shows the additional flat background, (c) One-dimensional energy cut obtained by integrating intensity in the range of 


1.2 A ^ < Q < 1.4 A The energy resolution is shown I 
integrating intensity in the range of 3.5 meV < E < 4.5 meV. 


(Cu-CHD), meaning that the origin of the excitation is 
a cluster of O 2 molecules. 

The excitation at ?ia;~2.4 meV exhibits stronger in¬ 
tensity at high Q. This means that its origin is nuclear 
lattice or cluster rather than magnetic one. We will ex¬ 
clude the excitation in the analysis section. Similarly we 
consider that the excitation at fia;~2.4 meV in O. 3 O 2 - 
(Cu-CHD) would have nuclear origin. 


IV. ANALYSES 

Combination of synchrotron x-ray diffraction and 
maximum entropy method/Rietveld analysis revealed 
that the adsorbed O 2 molecules form dimers in the 
nanopore. [ 2 ^ We, therefore, start our analysis on the 
magnetic excitation at i?~4 meV in 0.3O2-(Cu-CHD) 
from the 5 = 1 spin-dimer model shown in Fig. [SKa). We 
assume that the molecular axes of two O 2 molecules are 
parallel. The distance between the oxygen molecules, d, 
is fixed with the reported value 3.22 A.[^ The magnetic 


horizontal bar. (d) One-dimensional Q cut obtained by 


anisotropy axis is along the molecular axis, which is 
defined as 2 ; - axis. The spin Hamiltonian is 

nd^mer = J^i • ^2 + ^ (5f )^ (1) 

I 

The neutron cross section within the Born approx¬ 
imation [ 2 ^ is calculated by the diagonalization of 
the Hamiltonian. The powder average was, then, 
performed to compare the calculation with the ex¬ 
periment. We used the magnetic form factor 
^iQ) = drd9r^sin^9e~'^’^^^[cosh{2bRoCOs9) — 

2^jgiQrCOS(e-/3) described in Ref. HO- Here, we used the 
internuclear distance Rq=\.21 A and a constant 6=4.1 
A“^. /3 is the angle between the molecular axis and the 
scattering vector. 

Since the calculated neutron cross section of a spin 
cluster is a dispersionless delta function, the peak width 
of the excitation should be resolution limited. In con¬ 
trast the excitation, which is guided by the green curve 
in Fig. [1] (e), exhibits the larger width than that of the 
experimental resolution. The shape of the excitation is 
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(a) Bare Cu-CHD raw data (b) 2.0O2-(Cu-CHD) raw data 

£,=3.135 meV £,=3.135 meV 



^ ^ ^ I —u \j.yj I I I I 

0.0 0.5 1.0, 1.5 2.0 0.0 0.5 1.0, 1.5 2.0 

Q (A-i) Q (A-i) 
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I 
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Q (A-') 



(c) 




FIG. 3: (color online) (a,b) Contour map of the raw neutron spectrum of (a) bare Cu-CHD and (b) 2.0O2-(Cu-CHD) by 
using ifi=3.135 meV. (c) One-dimensional energy cut in the low-energy region obtained by integrating intensity in the range 
of 1.0 A ^ < <3 < 1.2 A ^ for raw 2.0O2-(Cu-CHD) and bare Cu-CHD by using i5i=3.135 meV. (d,e) Contour map of the raw 
neutron spectrum of (d) bare Cu-CHD and (e) 2.0O2-(Cu-CHD) by using Ei=7.743 meV. (f) One-dimensional energy cut in 
the high-energy region obtained by integrating intensity in the range of 1.2 A < Q < 1.4 A ^ for raw 2 . 0 O 2 -(Cu-CHD) and 
bare Cu-CHD by using £'^=7.743 meV. 


rather symmetric, and the asymmetry due to Van-Hove 
singularity is not observed. This suggests that the finite 
width is ascribed not to the powder-averaged dispersion 
due to an interdimer coupling but to the distribution of 
the intradimer exchange constants of isolated dimers. We 
consider that multiple minima for O 2 molecules are there 
inside the nanopore of Cu-CHD, leading to randomness 
of the configuration of O 2 dimer. Then the spin system 
is considered as a group of dimers having different Js. 
We assume that the distribution of the Js is Gaussian 
function. The free parameters for the fitting to the data 
are, thus, the mean value of the exchange interaction 
redefined as J, and the standard deviation of the distri¬ 
bution a. In this analysis, the uniaxial anisotropy D is 
not very important because the magnitude of a is much 
larger than that of D. We fixed D = 0.41 meV that is 
independently determined in the analysis of 2.0O2-(Cu- 
CHD) which will be explained later. We also consider the 
additional background with the form of A{huj)Q^ + B. 

For the fitting we use the data of one-dimensional 
energy cut in Fig. He) and one-dimensional Q cut in 
Fig. H^)- In the first step we fit the one-dimensional 
energy cut by three Gaussians to obtain the peak energy 
and the energy width. We introduce the additional back¬ 
ground with the phonon-like form i.e., A„exp[—((So; — 


TABLE I: Parameteres obtained for the 0.3O2-(Cu-CHD) 
sample from INS experiment (upper) and those for the 
0.22O2 -(Cu-CHD) sample from magnetization measurement 
(lower). 


J (meV) d (A) g (meV) D (meV) N{%) AE 


INS 

4.15 

3.22 2.16 

0.41 

- 

- 

M Dimer 

4.15 

1.77 

0.41 

90 

1.15J 

M Monomer 

- 

- 

0.41 

10 

- 


fujJn)/V^^n)^]Q^ {n = 1,2,3) and the additional flat 
background B. Here hiOn are the peak energies ob¬ 
tained in the first step. In the second step we fit the 
sum of the neutron cross section of the group of dimers 
and the additional background, ^,^(A„exp[—((lio; — 
hiJn)/+ B to the data of one-dimensional 
Q cut and one-dimensional energy cut. Fitting parame¬ 
ters are J, a, scale factor for the cross section of group 
dimers, A„(n = 1, 2, 3), and B. 

The calculated cross sections with the parameters sum¬ 
marized in Table I are shown by by the green solid curves 
in Figs. He) and[2Dd). The blue curves are the additional 
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(a) 2.0O2 -(Cu-CHD) subtracted data (b) 2.0O2-(Cu-CHD) subtracted data 



(C) 




FIG. 4: (color online) Neutron intensity for 2.0O2-(Cu-CHD) after the background subtraction. The green curve is the magnetic 
excitation, the blue curve is the additional background having Q dependence, and the red curve is the sum of green and blue 
curves. The dashed line shows the background, (a) One-dimensional energy cut in the low-energy region obtained by integrating 
intensity in the range of 1.0 A < Q < 1.2 A . The energy resolution is shown by the horizontal bar. (b) One-dimensional 
energy cut in the high-energy region obtained by integrating intensity in the range of 1.2 A ^ < Q < 1.4 A (c) One¬ 
dimensional Q cut obtained by integrating intensity in the range of 0.3 meV < E < 0.5 meV. The energy resolution is shown 
by the horizontal bar. (d) One-dimensional Q cut obtained by integrating intensity in the range of 3.5 meV < E < 4.5 meV. 


background, and the red solid curves are the sum of green 
and blue curves. The data are reasonably reproduced 
by the calculation, meaning that the magnetic excita¬ 
tion of the 0.3O2 -(Cu-CHD) at T = 4.5 K is explained 
by the spin dimers with normally distributed exchange 
constants. 

The eigenenergies of a single dimer calculated by using 
the parameters of J and D in Table I are shown by the 
black lines in Fig. EJc). The ground state, Eg, is non¬ 
magnetic singlet. The first group of the excited states, 
Eist, is composed of nearly triplet states with S'=l, which 
are lifted by the single-ion anisostopy D. The second 
group, E 2 nd, is composed of nearly quintet states with 
S=2. The transition between Eg and Ei^t is probed as 
the excitation at lEj = 4 meV in the INS spectrum. 

In the INS spectrum in 2.0O2-(Cu-CHD) we observed 
a couple of magnetic excitations at i?~0.4 meV and 4 
meV. The latter excitation looks similar to that observed 
in 0.3O2 -(Cu-CHD), but the simple dimer model does not 
explain the 0.4 meV excitation. In fact the triplet state 
is the ground state in case that the cluster is composed 
of the odd number of Heisenberg spins, and, therefore, 
the transition between the triplet states lifted by the 
single-ion anisotropy D can be probed as the low en¬ 


ergy excitation. In the monomer case, the cross section 
decreases monotonically with increasing Q because it de¬ 
pends solely on the magnetic form factor. In Fig. [3] (g) 
the cross section exhibits a peak at Q~1.2 A“^, meaning 
that the cluster is composed of the multiple spins. We, 
thus, consider the trimer of O 2 molecules shown in Fig. 
[5Kb). We assume that the molecular axes are parallel. 
The trimer has two different inter-molecular lengths, di 
and c? 2 , and correspondingly two different exchange con¬ 
stants, Ji and J 2 - The spin Hamiltonian of a trimer is, 
thus, 

'^trimer — Ji 5 i • S2 + J2S2 -Ss + dY, (sn^ ■ ( 2 ) 

I 

The neutron cross section is calculated by the diagonal- 
ization of the Hamiltonian. The free parameters are the 
exchange interactions Ji, J 2 , the uniaxial anisotropy D, 
and one of the inter-molecular distances c? 2 - Another 
distance, c?i, is fixed to 3.22 A. [ 2 ^ We assume that the 
width of the magnetic excitation at A~4 meV is due to 
the normal distribution of the exchange constant Ji. We 
define that the standard deviation cri of the distribution, 
and redefine the mean value of the exchange constants as 
Ji. Both of them are the free parameters for the calcula- 
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(e) H configuration (f) X configuration (g) S configuration 



The energy level of a single trimer calculated by the 
obtained parameters and the mean exchange constant Ji 
are shown by the black lines in Fig. [5td). The eigen¬ 
states of the trimer model are separated into the three 
groups. The group having the low energy, Eq, is 5" = 
1 states lifted by the single-ion anisotropy D. They are 
composed of 3 states and the ground state is Sz=0. The 
first group of the excited states, is composed of 9 
states, and the second group, i? 2 nd, is composed of 15 
states. The excitation at i?~0.4 meV corresponds to the 
transitions among the Eq states. The excitation at i?~4 
meV corresponds to the transitions between Eq states 
and Eist states. 

The Q dependence of the excitation at 4 meV in 2 .OO 2 - 
(Cu-CHD) is similar to that in 0 . 3 O 2 -(Cu-CHD), since 
they are physically quite similar. If J 2 is zero and the 
system is composed of isolated dimers and monomers, the 
transition between Eq and Eigt in the system is exactly 
the same as that in dimers system. It is, therefore, rather 
difficult to exclude the possibility that some amount of 
dimers are included in 2 . 0 O 2 -(Cu-CHD). Even though 
the intensity comparison of the excitation at 0.4 meV, 
which is purely from spin timer, and that at 4 meV tells 
us the ratio of dimers to trimers in principle, the sepa¬ 
rated data set in the different experimental setups makes 
the comparison difficult. 


FIG. 5: (Color online) (a) O 2 configuration used in the cal¬ 
culation of the dimer model, (b) O 2 configuration used in the 
calculation of trimer model, (c) Black lines show energy lev¬ 
els calculated by conventional spin Hamiltonian of 5=1 dimer 
model. Red lines show the reduced energy levels to repro¬ 
duce the magnetization, (d) Black lines show energy levels 
calculated by conventional spin Hamiltonian of 5=1 trimer 
model. Red lines show the reduced energy levels to reproduce 
the magnetization, (e) O 2 dimer in H-configuration. (f) O 2 
dimer in X-configuration. (g) O 2 dimer in S-configuration. 


tion. The standard deviation (T 2 of J 2 is fixed to 0. The 
additional background is also considered in a similar way 
to the analysis in 0 . 3 O 2 -(Cu-CHD). 

The calculated cross sections are shown by the solid 
curves in Figs. He)®, and the obtained parameters 
are shown in Table II. The green curve is the magnetic ex¬ 
citation, the blue curve is additional background, and the 
red curve is the sum of green and blue curves. The data 
is reasonably reproduced by the spin trimers of which 
the main exchange constants are normally distributed. 
We note that the ratio of the peak intensity at Q~1.2 
to the amplitude of the intensity modulation in 
Fig- El (g) is tuned by the ratio of Ji to ^ 2 - The value 
of d 2 is larger than that of di, and the value of J 2 is 
much smaller than that of Ji. This means that the O 2 
trimers in the high oxygen concentration sample 2 .OO 2 - 
(Cu-CHD) are regarded as the O 2 dimers weakly coupled 
to O 2 monomers. 


V. DISCUSSION 

The magnetic held dependence of the energy level of 
Heisenberg 5=1 spin dimer is shown in the upper panel 
of Fig. [6] (b), and the corresponding magnetization curve 
is shown in the lower panel. The magnetization curve ex¬ 
hibits 1/2 plateau at He 5 77 < 2He with He = J/{gII b): 
where the ground state is |5, 5^) = |I, 1). In contrast 
in 0.22O2 -(Cu-CHD) the plateau is absent [l^ as shown 
in Fig. [71(a), and it seems more likely two-level system. 
Anomalous energy level of 5 = I spin dimers in the ad¬ 
sorbed oxygen system was initially discussed in the tem¬ 
perature dependence of the neutron intensity in a differ¬ 
ent metal complex, CPL-I.[I^ The intensity decreased 
more rapidly with the increase of the temperature com¬ 
pared with conventional 5 = 1 spin dimer model, sug¬ 
gesting the lowering of the energy level of the 5 = 2 
quintet state and the proximity to the triplet state of 
5=1. The energy level was explained based on the 
calculation including inter-molecular potential of 02 . 0 ] 
For the dimer of O 2 molecules, both the ground state 
with 5=0 and the hrst excited states with 5=1 have an 
H-configuration shown in Fig. HKe) in the equilibrium 
molecular positions. On the other hand, the second ex¬ 
cited state with 5 = 2 has an X-configuration shown in 
Fig. 0(1) of which the energy is lower than that in the 
H-conhguration. The energy difference between 5 = 0 
and 5 = 1 in both H configuration and that between 5 
= 1 in H configuration and 5 = 2 in X conhguration is 
almost the same, resulting in the anomalous energy level. 















































TABLE II: Parameteres obtained for the 2.0O2-(Cu-CHD) sample from INS experiment (upper) and those for the 1.1102-(Cu- 
CHD) sample from magnetization measurement (lower). 


Ji (meV) di (A) CTi (meV) J 2 (meV) ^2 (A) cj 2 (meV) D (meV) N{%) AE 


INS 3.94 

3.22 1.54 

0.62 

3.5 

0 

0.41 

- 

- 

M Trimer 3.94 

2.52 

0.62 

- 

0 

0.41 

33 

1.3J 

M Dimer 4.15 

3.54 

- 

- 

- 

0.41 

67 

1.3J 


(a) Monomer (b) Conventional dimer (c) Modified dimer 




FIG. 6: The field dependence of the energy levels and corresponding magnetization curve in S' = 1 spin cluster systems, (a) 
Monomer spin with a single-ion anisotropy D. (b) Conventional Heisenberg spin dimer, (c) Heisenberg spin dimer with the 
reduced energy levels of S = 2 states, (d) Isolated monomer and dimer, (e) Conventional trimer where a dimer is weakly 
coupled to a monomer. Single-ion anisotropy is also considered, (f) The trimer with the reduced energy levels of higher energy 
states. 


In the magnetization and x-ray study in O. 22 O 2 - 
(Cu-CHD), the lowering of the S = 2 state was ex¬ 
plained in the similar manner. According to their Ri- 
etveld/MEM analysis S'-geometry as shown in Fig. Elg) 
was realized for S' = 2 state. It was estimated that 
the first energy-gap between S = 0 and S = I is 3.7 
meV and the second one between S = I and S = 2 


was 4.7 meV. The distribution of the exchange constant 
and the single-ion anisotropy term were not considered 
there. In the present paper we consider the magneti¬ 
zation data of both 02 /(Cu atom)=0.22 and 1.11 re¬ 
ported by W. Mori et al. [3 Our discussion is based 
on the spin monomer/dimer/trimer model with the nor¬ 
mally distributed exchange constant and the single-ion 
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anisotropy obtained from the INS spectra. We also as¬ 
sume that the eigenenergies of higher states calculated 
by conventional spin Hamiltonian are reduced because 
of the spin-dependent molecular potential in the oxygen 
dimer. [GHoI 


We consider the following Hamiltonians 


"^mono 

1 1 

(3) 

^dimer 

= Js^-s^ + DYisn^-gt^BHY 

1 1 

Sf{4) 

^trimer 

= Ji5i • 52 + J 2 S 2 -Ss + dY )" 

1 



- gfioHYsf 

(5) 


i 


where the direction of the magnetic field, H, is defined as 
z'-axis. For all models, we fixed D=QA\ meV. The mag¬ 
netizations were calculated using the above Hamiltonians 
having the distributions in J and Ji. The powder average 
was performed. 

The reported magnetization in 0 . 22 O 2 -(Cu-CHD) ex¬ 
hibited a small increase below 5 T and a large increase 
at around 20 T with increasing H . The former increase 
is not explained by S' = 1 dimer model but by S = 1 
monomer model. As schematized in Fig. IHIa) the mag¬ 
netization of the monomer at T = 0 K exhibits a step 
function of which the jump held is Ho = D/{gfiB)- Use 
the value of D = 0.41 meV and Ho = 3.5 T is obtained. 
At hnite temperature the jump is smeared. The held 
scale of the initial increase of the magnetization is, thus, 
consistent with that of the single-ion anisotropy D = 0.41 
meV. We, therefore, assume that S = 1 oxygen molecule 
monomers were trapped in some local minima in the 
nanopore or surface of the grain of the crystal in 0 . 2202 - 
(Cu-CHD). The ratio of O 2 molecules in the monomer, 
N, is estimated to be 10 % of the total O 2 molecules from 
the magnetization plateau at 8 T < H < 15T. The small 
number of the monomers would be the reason why the 
excitation was not observed in INS spectrum. 

The red lines in Fig. [5jc) shows the energy levels of the 
spin dimer of which the energy of the quintet S = 2 state 
is lowered compared with the conventional spin dimer. 
The field dependences of the conventional and modified 
states are shown in the upper panel in Fig. [Sib) and (c). 
The single-ion anisotropy is omitted for the simplicity 
here. In the modified scheme the energy of the = 
2 crosses at the lower field than that of = 1 does, 
meaning that the magnetization curve is a single-step 
function as shown in the lower panel. By introducing the 
distribution of the exchange constant, the critical field is 
smeared and the experimentally obtained magnetization- 
curve is reproduced. We assume that the inflection point 
of the experimental magnetization curve is the critical 
field of the spin dimer having the mean value of J, 4.15 
meV, in Table I. The field of the inflection is estimated 
to be 34 T, leading to that the decrease of the energy 
of the S = 2 state, AE, is 1.1 J. Considering the single- 


(a) (b) 



«(T) H(J) 


FIG. 7: (Color online) Comparison between the magne¬ 

tization reported by W. Mori et ab[I^ and the calculated 
magnetization. The black symbols show the reported data, 
and the blue and red solid lines represent the calculation (see 
text), (a) Magnetic field dependence of the magnetization 
for Cu-CHD adsorbing less oxygen molecules, (b) Tempera¬ 
ture dependence of the magnetic susceptibility for Cu-CHD 
adsorbing less oxygen molecules, (c) Magnetic field depen¬ 
dence of the magnetization for Cu-CHD adsorbing more oxy¬ 
gen molecules, (d) Temperature dependence of the magnetic 
susceptibility for Cu-CHD adsorbing more oxygen molecules. 


ion anisotropy we found that 1.15 J is more appropriate 
value. The distribution of the J is manually changed 
to find the best consistency. It is find that 1.77 meV 
is the best standard deviation of J. The magnetization 
curve calculated by using the parameters in Table I is 
indicated by the red solid curve in Fig. [ZKa). The data is 
reasonably reproduced by the calculation. The blue solid 
curve is the magnetization for cr = 0, meaning that the 
distribution of J is required. The component of dimer 
and that of monomer is indicated by green solid curve 
and black dashed curve in Fig. [T^b). 

The lowering of the S = 2 state is not probed in the 
INS profile because the transition between 5 = 0 and S = 
2 is forbidden. The transition between 5=1 and 5 = 2 
is not forbidden but is not observed because of the small 
probability of 5 = 1 state at the low temperature. We 
note that the excitation at E~2A meV observed in INS 
is not probed in the magnetization. The result supports 
our consideration that it is not from the magnetic origin. 

The magnetization of 1.1102-(Cu-CHD) measured by 
Mori et al. [l^ is indicated by solid symbols in Fig. 
[TKc). To understand the qualitative behavior of the 
5=1 trimer composed of a dimer weakly coupled to 
a monomer, it is instructive to consider an isolated 
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monomer and dimer. The energy level and the corre¬ 
sponding magnetization curve is shown in Fig. [Hd). The 
magnetization exhibits three-step staircase where the 
critical fields are Hi^c and 2Hi c with Hi c = Ji/{g^B)- 
The ground state is S' = 1 triplet and the magnetization 
has 1 (//xs/trimer even at 0 T. By introducing the weak 
coupling between dimer and monomer and the single¬ 
ion anisotropy D, the degenerated states are lifted and 
the energy levels and the corresponding magnetization 
curve are modihed as shown in Fig. [TJ^e). Here the high¬ 
est states are |S, S^) = |3, ±3) and the energy is about 
3Ji. In order to make the magnetization a two-step stair¬ 
case, it is necessary to lower the energies of the group of 
the higher states at about 3Ji by Ai? as shown in the 
upper panel in Fig. El^f). By introducing the distribution 
of Ji the step function is smeared. We tried to reproduce 
the experimental data by the trimers model but we failed. 
Instead the data was reproduced by the sum of the dis¬ 
tributed dimers and the distributed trimers as indicated 
by red solid curve in Fig. [7{c). The ratio of O 2 in dimers 
is 67% of the total adsorbed O 2 . In Fig. [TJd) the magne¬ 
tization of the dimers and trimers are indicated by green 
and black solid curves, respectively. The parameters used 
in the calculation are summarized in Table II. 

The concentration of O 2 in the sample used in magne¬ 
tization measurement is l.ll[mole/f.u.], and it is smaller 
than that in INS measurement 2.0[mole/f.u.]. The less 
amount of O 2 would be the reason for the substantial 
number of dimers in the magnetization data. The data, 
however, indicates that some amount of dimers may be 
included also in the sample of 2.0O2-(Cu-CHD) used in 
the INS experiment. The existence of several potential 
minima in the nanopore may be a possible reason for the 
inclusion of both dimers and trimers. 


VI. CONCLUSION 


We performed INS measurements on Cu-CHD adsorb¬ 
ing O 2 molecules with the low and high concentrations 
to identify the magnetism of the 02 -based magnet re¬ 
alized in the nanopores. The INS spectra of the two 
different samples are explained by different spin Hamil¬ 
tonians, spin dimers and spin trimers, meaning that the 
spin system can be controlled by the concentration of 
O 2 molecule in the 02-based magnet. It is found that 
the magnetic excitation is broad and, therefore, the sys¬ 
tem is composed of a group of clusters having normally 
distributed exchange constants. This indicates that the 
supercrystal of the oxygen molecules is not perfect. By 
using the parameters obtained in INS and by assum¬ 
ing the reduction of the higher energy states due to the 
non-negligible spin-dependence in the molecular poten¬ 
tial, magnetization curves are explained in quantitative 
level. 
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